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WORK  PROGRESS 

Progress  during  the  second  quarter  is  presented  in  this 
report  by  principal  tasks.  Areas  of  current  investigation  and  plans 
for  the  third  quarter  are  also  presented. 
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Using  an  ensemble  of  seismic  array  network  data 
to  be  furnished  by  AFTAC,  investigate  bodywave 
noise  on  a  coherent  worldwide  basis.  Investigate 
interarray  equalization  problems.  Study  methods 
of  combining  the  subarray  output  for  network 
signal  extraction.  Investigate  the  capabilities  of 
a  worldwide  network  for  resolving  events  closely 
spaced  in  time  and  space. 

Data  required  for  the  network  signal  and  noise  studies  under 
Task  A  are  routinely  received  from  SDL.  The  data  are  digitised  from 
FM  magnetic  tape  recor.  lings,  demultiplexed  and  reformatted  for  input 
to  IBM  equipment  by  SDL.  DAC  playbacks  of  the  digitized  records 
are  also  being  provided  by  SDL.  A  new  IBM  7044  program  has  been 
written  and  checked  out  which  provides  translation  of  the  data  into 
acceptable  format  for  use  with  existing  7044  software.  All  digital  tapes 

received  from  SDL  to  date  have  been  satisfactorily  translated  with  this 
progra  i\. 

Data  for  six  events  and  five  noise  samples  have  been 
received  from  SDL,  with  part  of  the  ten  requested  stations  for  each 
event.  The  remaining  nine  events  have  been  requested;  and  in  addition, 
long-period  data  sampled  at  a  K- sample /sec  rate  have  been  requested 
for  these  events.  Appropriate  long-period  and  short-period  calibrations 
have  also  been  requested.  Instrument  response  curve  and  array  con¬ 
figuration  information  have  been  obtained  and  when  further  information 
pertaining  to  SDL  filter  response  curves  is  received,  all  noise  spectra 
will  be  converted  to  units  of  absolute  ground  motion.  Meanwhile, 
analysis  of  calibration  data  already  received  is  underway. 
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Long-period  data  have  been  requested  for  two  noise 
samples  for  which  short-period  data  are  already  available.  These  data 
will  be  analyzed  to  determine  whether  long-period  or  short-period 
data  are  more  suitable  for  noise  characterization  (i.e.,  study  of  high- 
velocity  noise)  at  the  network  level.  Additional  noise  sample  requests 
are  withheld  pending  the  results  of  this  evaluation. 

An  initial  analysis  of  the  noise  data  has  been  conducted  to 
establish  the  seismic  validity  of  the  recorded  data.  Power  density 
spectra  of  five  noise  samples  have  been  computed  for  one  channel  each 
of  all  available  stations.  These  spectra  indicate  seismically  valid 
data  up  to  2.0  cps,  with  seismic  noise  20  to  40  db  above  instrument 
noise  in  this  band.  Above  2.5  cps,  most  spectra  exhibit  narrow  peaks 
probably  related  to  instrument  noise. 

The  various  array  configurations  have  been  examined,  and 
spectral  windows  and  straight  sum  response  for  each  array  have  been 
computed.  A  comparison  of  several  methods  for  computing  f-k  spectra 
is  underway  with  the  most  suitable  method  subsequently  used  in  the 
high  velocity  noise  analysis. 

Several  techniques  for  velocity  filtering  are  under  considera¬ 
tion,  and  a  comparison  is  bein*.  made  of  the  relative  effectiveness  for 
signal  extraction  of  straight  sum,  measured  noise  Wiener  MCF  and 
adaptive  maximum  likelihood  filtering  techniques.  Primary  emphasis 
is  being  placed  on  comparing  the  results  of  the  adaptive  approach  to 
the  standard  methods  in  terms  of  performance  and  cost. 

Although  explosion  data  are  required  for  an  optimum  investiga¬ 
tion  of  interstation  signal  equalization,  a  preliminary  investigation  of 
equalization  problems  at  the  network  level  has  begun  using  earthquake 
data  only. 


Task  B 


5  July  1967 
Page  4 


Continue  investigations  of  multi-element  system  studies 
to  determine  possible  new  combinations  o';  senstfrs  for 
noise  reduction.  Study  methods  of  specifying,  for 
given  noise  fields,  the  optimum  multisensor  system 
for  noise  reduction  with  the  desired  result  a  set  of 
guidelines  for  array  design.  These  guidelines  should 
include,  but  not  be  limited  to,  the  type  of  sensors 
required  and,  for  an  array,  the  s?ze  and  geometry, 
subject  to  the  constraints  of  practicability. 

1.  Program  for  Computing  Theoretical  Crosspower  Spectra  and/or 
Crosscorrelation  Functions 

Calculation  of  crosspower  spectra  from  the  solid  cone 
model  when  horizontal  seismometers  are  present  has  been  found  to  be 
impractical  in  the  framework  of  this  program.  It  is  possible  to  handle 
horizontal-vertical  arrays  for  a  model  which  appears  in  f-k  space  as 
the  surface  of  a  cone.  The  solid  cone  model  can  then  be  approximated 
by  a  sum  of  a  group  of  "surface-of-a-cone"  models.  This  surface-of- 
a-cone  model  is  to  be  distinguished  from  a  degenerate  case  of  surface 
wave  propogation  which  assumes  the  same  geometric  shape  in  f-k 
space.  The  first  arises  from  compressional  waves  emerging  from  depth 
m  a  homogeneous  half  space.  The  waves  are  assumed  to  approach  the 

surface  with  a  single  angle  of  incidence  and  are  equally  likely  to  come 
from  any  direction.  The  second  model  represents  non-dispersive 
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Rayleigh  wave  propogation  which  can  exist  in  a  homogeneous  half- 
space.  The  second  model  as  well  as  dispersive  surface  wave  modes 
can  be  handled  in  the  program. 

Actual  programming  has  been  delayed  pending  completion 
of  the  IBM  360  installation.  At  this  time,  it  appears  that  programming 
can  begin  in  one  month.  Current  effort  is  devoted  to  compressional 
wave  models  which  arise  in  the  presence  of  a  layered  media. 

2.  WMSO  Noise  Field  Study 

Analysis  of  the  1962  noise  sample  from  WMSO  has  been 
completed.  A  special  report  describing  the  results  of  this  study  is 
in  preparation  and  will  be  submitted  during  the  month  of  July. 

3.  Augumentation  of  the  WMSO  Array  and  Horizontal  Sensors 

Based  on  the  noise  study  results,  a  recommendation  has 
been  submitted  for  the  emplacement  of  horizontal  seismometers  in 
the  array  at  WMSO.  This  would  permit  experimental  study  of  such 
sensors'  value  in  the  enhancement  of  body  wave  detection  in  the  presence 
of  surface  wave  noise.  The  relatively  high  level  of  surface  wave 
energy  at  WMSO  makes  this  a  desirable  location  for  the  study. 

Task  C 

Theoretically  investigate  methods  of  implementing  contin¬ 
uously  adaptive  systems  for  application  to  time -varying 
noise  fields  and  postdetection  processing.  Any 
system  that  can  be  simulated  off-line  should  be 
evaluated  using  suitably  characteristic  data. 

One  way  of  viewing  and  perhaps  developing  the  adaptive 
algorithm  is  as  a  classical  iterative  matrix  inversion  technique.  The 
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classical  technique  is  modified  by  letting  the  matrix  to  be  inverted, 
along  with  the  solution,  vary  with  time  or  iteration. 

In  solving  for  x  in  the  system  of  equations  Ax  +  r  =  0, 
one  may  start  from  an  approximate  solution  x  =  Ufc  and  proceed  in 

the  direction  r  =  A  +  r  a  distance  -  \r  so  that 

uf  uf  u. 


u 


t+1 


For  sufficiently  small  X,  X  <  2/X^  ,  the  sequence  u(  will  converge 
to  the  solution  x. 

An  autoregressive  process  of  degree  p  can  be  represented 
by 


X 


t+1 


=  a 


X 


(P) 


+  e. 


t  +  1 


where 


(X  .  X 


t’  ~t-l’ 


•Xt-P+1> 


and  a  is  a  vector  of  fixed  coefficients  and  the  are  independent 

N(0,  1)  variables.  It  is  known  that  a  satisfies 
Ra  =  y 


where 
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Y  =  E  Xp}  X 

t  t+1 


At  any  point  in  time, 


'  E  reasonable  estimate  of  R,  R^  is  given  by 

Rt  =  MRt-l  +  (1"^  X|P)  (X(P^)T 

^  t 


and,  similarly, 


Yt  =  UYt.l  +  -M)  xjp)  Xt+]  . 

The  update  procedure  for  ot  is 

“t  =  Vi  +  UYt-Etat.i> 

and  if  we  let  u ,  0,  then  the  adaptive  procedure 
follows. 


we  have  been  using 


One  reason  for  looking  at  the  adaptive  algorithm  from  this 
pornt  of  view  is  that  the  same  reasoning  can  now  be  applied  to  general 
signal  extraction  problems  with  arbitrary  but  unknown  signal  models. 
The  general  form  of  the  equation  to  be  solved  is 

(Ro  +  Bn)  a  =  ys  . 

The  estimates  of  ^  and  a  are  again  indexed  with  t  to  get  the  adaptive 
algorithm. 


°t+i  =at  - 


d  =  [Bs  +  «n  >t  ]  at 
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Ihere  are  numerous  better  ways  to  implement  the  above  equations 
through,  approximations  of  and  these  are  under  investigation  at 
present.  Also,  it  is  clear  that  each  of  the  classical  matrix  inversion 
iterative  techniques  corresponds  to  a  potential  adaptive  processing 
algorithm. 

Time  Constant  Considerations 

It  is  desirable  to  be  able  to  relate  the  parameter  \  to 

the  effective  amount  of  data  that  is  being  used  to  design  the  set  of 

*  ^ 
filter  coefficients  at  each  point.  Defining  the  sequence  Ct  by 

*  '  (P) 

at  =  at +  \  et  x. 


where 


\  =  11  [(X[P)  )T  (X|P)^  3 


£  =  X  -  aT  x<P> 

et  Xt+1  at  Xt 


and  a^_  is  the  adaptive  algorithm  sequence  of  filter  vectors,  the  sequence 

* 

ttt  *S  'l  ,  ^a-*-8e  gain"  set  of  coefficients  since 


❖ 

e  =  X. 


*  T  (P) 

at  xi  ■  "• 


However,  the  important  thing  to  consider  at  present  is  that  e  =  0 
implies  that  the  time  constant  of  the  with  respect  to  the  data  is 

very  short.  Now  consider  the  sequence 


00 


V  -i/T 
L  e 
i=l 
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8t_!  +  (1-e  '1/T)  a* 


t-1 


_  J. 

=  Pt-1  +  (1-e  )  (a*  J  -  tt  i) 


and  it  follows  that  3(  is  the  a,  defined  by  the  adapts  algorithm 
\  =  (1-e  )  X  =  c\  . 

The  time  constant  of  the  sequence  <xt  with  respect  to  the  data  is 

therefore  taken  to  be  t  so  that  an  "equivalent"  boxcar  section  of  data 
of  length  N  is  given  by 

N  -  (2-c)/c. 

Task  D 


Investigate  the  effects  and  methods  of  reducing  locally 
generated  noise.  The  effects  of  such  non-plane  wave 
fields  on  multichannel  filter  design  should  be  evaluated. 

Adaptive  processing  of  four  sections  of  seismic  multi¬ 
channel  data  has  been  completed  and  a  report  of  the  results  is  in 
prcgress.  At  present  it  appears  that  adaptive  prediction  and 
adaptive  signal  extraction  by  maximum- likelihood  filtering  can  be 
performed  to  almost  duplicate  the  Wiener  results  for  actual  seismic 
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data.  In  addition  to  being  logically  simpler  to  implement  and  less 
expensive,  the  adaptive  approach  appears  to  be  effective  in  pro¬ 
cessing  time -varying  data. 

No  further  work  on  near-array  noise  source  problems 
has  been  done.  Work  was  begun  in  June  that  makes  use  of  programs, 
data,  and  experience  accumulated  in  the  work  performed  under 

Task  C. 

Task  E 


Continue  studies  of  the  instrumental  equalization 
problem.  Apply  any  new  techniques  available  for 
studying  instrumental  equalization  and  evaluate 
the  effectiveness  of  such  techniques. 

1.  Statistical  Noise  Prediction 

Signal  extraction  for  a  p-channel  array  frequently  reduces 

to  a  noise  estimation  problem  (maximum  likelihood  processing,  for 
example)  in  which  the  noise  crosscorrelations  are  estimated  to  determine 
an  estimated- filter  which  could  be  used,  for  instance,  to  predict  the 
noise  on  one  channel  f'om  the  noise  on  the  other  p-1  channels.  The 
performance  (mean-square  prediction  error)  of  such  estimated-filters 
has  been  compared  with  the  optimal  Wiener  filter  error.  The  ratio 
of  the  estimated-filter's  performance  to  the  optimal  filter  performance 
is  independent  of  the  true  channel  crosscorrelations.  Thus,  the  degree 
of  channel  correlation  does  not  affect  the  relative  performance  of  the 

estimated-filter. 
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Similarly,  the  ratio  Of  the  estimated-filter's  performance 
on  the  design  points  to  the  optimal  Wiener  error  is  independent  of 
the  channel  crosscorrelations.  Thus,  the  relative  false  gain  fo»*  a 
highly  correlated  array  is  no  worse  than  for  a  less  correlated  array. 

If  the  noise  is  assumed  Gaussian,  the  statistics  of  these 
invariant  ratios  can  be  determined.  Such  calculations  tell  the  designer 
how  well  the  estimated  filter  can  be  expected  to  perform.  For 
example,  a  12-channel  filter  based  on  90  independent  samples  will 
perform  within  1  db  of  the  optimal  (but  unknown)  Wiener  filter  with 
99  percent  confidence.  Such  computations  are  extremely  useful  in 
determining  how  many  samples  are  needed  to  achieve  some  desired 
relative  performance  or  in  evaluating  the  performance  to  be  expected 
from  an  existing  system. 

A  special  report  which  derives  and  explains  these  results 
is  being  prepared.  Future  effort  will  be  devoted  to  extending  these 
results  to  other  related  multichannel  array  problems. 

2.  Long-Noise  Sample  Study 

Thirteen  noise  samples,  from  1  to  4  hours  long,  have  been 
selected  initially  from  LASA  Develocorder  film  recorded  in  December 
1966  and  April  1967.  Playbacks  of  two  of  the  better  samples  are 
being  requested  to  satisfy  data  requirements.  A  five  to  six  hour 
noise  sample  is  still  desired  and  will  be  obtained  if  it  becomes  avail¬ 
able. 
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